Small-angle neutron scattering by spatially inhomogeneous ferromagnets with a nonzero average uniaxial anisotropy

Macroscopic spin-misalignment small-angle neutron scattering cross sections and response functions are computed analytically for a material with global uniaxial magnetic anisotropy (texture). The resulting expressions are tested against previously published experimental data.

Current micromagnetic SANS theory, however, deals only with macroscopically isotropic magnets (Honecker & Michels, 2013;Metlov & Michels, 2015;Michels et al., 2016). Such systems can have small local fluctuating magnetic anisotropies, but their magnitude and direction average out, so that the sample as a whole remains statistically isotropic. There is an abundance of such samples, but there are also plenty of material systems which do have a nonzero global anisotropy (texture) (Michels et al., 2017;Oba et al., 2020). Moreover, anisotropy can be induced in an originally statistically isotropic magnet artificially, e.g. by annealing it in the presence of a magnetic field, applying mechanical stress or subjecting it to severe plastic deformation (Herzer, 1997(Herzer, , 2013.
The purpose of this paper is to extend the micromagnetic SANS theory to systems with a global uniaxial anisotropy. Starting from the relevant micromagnetic Hamiltonian, we compute the spin-misalignment SANS cross sections analytically for many practically interesting cases of mutual orientations between the external magnetic field and the anisotropy axis. The results reveal new effects and quantitatively describe well known observations. We test the theory by applying it to previously published experimental SANS data.
The work is organized as follows: in Section 2 we introduce the nomenclature and the expressions for the SANS cross sections in terms of the Fourier image of the magnetization distribution of the sample. Section 3 contains the solution of the micromagnetic problem for the magnetization distribution of a sample with global anisotropy. When the magnetic field is perpendicular to the global anisotropy axis, there are two distinct cases: the small-field limit (when the direction of the average magnetization is determined by the balance of the external field and the anisotropy torques) and a simpler highfield limit (when the average magnetization is aligned strictly along the external field direction, but its small fluctuations are still influenced by the global anisotropy of the sample). After performing the averaging of the SANS cross sections over the orientations and realizations of the random magnetic parameter fluctuations in the scattering volume, we present the resulting macroscopic cross section expressions in Sections 4 and 5 in the low-and high-field limits for several chosen mutual orientations between the anisotropy axis and the external magnetic field. In Section 6 we apply our theory to the interpretation of unpolarized SANS data on a fieldannealed nanocrystalline bulk Vitroperm metallic glass and on a pure Ni sample subjected to severe plastic deformation. Finally, Section 7 summarizes our main results.

Magnetic small-angle neutron scattering
A scheme of a typical SANS experiment is shown in Fig. 1. An incident neutron from the incoming beam with a well defined energy and wavevector k is scattered by the sample and takes on a new wavevector k 0 . The difference q ¼ k 0 À k is called the scattering or momentum-transfer vector. Magnetic SANS experiments are usually performed in the presence of an applied magnetic field H. Out of all the possible mutual arrangements of the vectors k and H, two scattering geometries are usually employed: the perpendicular one (k ? H) and the parallel one (k k H). We can associate a Cartesian coordinate system with each of the geometries in such a way that their OZ axis coincides with the direction of the applied magnetic field H and either the OX [in the perpendicular geometry, Fig. 1(a)] or the OZ [in the parallel geometry, Fig.  1(b)] axis coincides with the wavevector of incident neutrons. In SANS, the vector q is assumed to lie in the detector plane and is usually parametrized as q ? ¼ qf0; sin ; cos g or q k ¼ qfcos ; sin ; 0g in the perpendicular and the parallel geometry, respectively.
Assume that the sample is a magnetically ordered substance (ferro-or ferrimagnet) at a sufficiently low temperature, far from the order-disorder phase transition (and the compensation temperature, if it is a ferrimagnet). Such a medium can be characterized by a nonzero local magnetization vector M, jMj ¼ M S being the saturation magnetization of the material. At the mesoscopic level, the spatial distribution of the local magnetization (magnetic texture) can be represented by a continuous vector field MðrÞ, where r is a spatial coordinate.
If the material is infinite, isotropic and uniform, its equilibrium magnetization will be oriented along the direction of the external field H. Otherwise, there will be some spatial variation of the magnetization, which will manifest itself, in particular, via the magnetic scattering of neutrons at nonzero scattering vectors q. The corresponding total (magnetic and nuclear) macroscopic scattering cross section dAE=d can be expressed as the sum dAE=d = dAE res =d þ dAE M =d, where the first term (residual scattering) is assumed to be magnetic field independent and the second term (spin-misalignment scattering) vanishes under the condition of magnetic saturation (very large external magnetic field magnitude). In the perpendicular (superscript ?) and parallel (k) scattering geometries the unpolarized spin-misalignment scattering cross sections can be expressed as follows (Michels, 2021): b H = 2.906Â10 8 A À1 m À1 is the magnetic scattering length per Bohr magneton, V denotes the scattering volume, the polar angles and of the scattering vector q in the detector plane are schematically depicted in the insets of Fig. 1, M ¼ fM X ; M Y ; M Z g is the magnetization vector inside the material with Cartesian components M X;Y;Z (in the coordinate system for a particular scattering geometry), lines above the symbols stand for the complex conjugate, and tildes denote the Fourier transform: The integral in (3) is taken over a representative cube of the material V ¼ L Â L Â L, corresponding to the coherence volume of the neutron beam, which is considered to be periodically repeating. The Cartesian components of q ¼ fq X ; q Y ; q Z g therefore take on all values which are integer multiples of 2=L. Because of this different definition of the Fourier transform, expressions (1) and (2) differ by the factor V 2 =ð8 3 Þ from the ones in the original work (Honecker & Michels, 2013). Thus, computing the SANS cross sections (1) and (2)

Micromagnetics of a weakly inhomogeneous yet globally anisotropic magnet
Micromagnetics (Brown, 1963) is based on the minimization of the total energy of the magnet E, whose local magnetic moments are subject to various well known magnetic interactions. Specifically, in this work, similarly to Metlov & Michels (2015), we assume that the total energy can be expressed as an integral over the energy density where CðrÞ is the position-dependent exchange stiffness, r r r ¼ f@=@X; @=@Y; @=@Zg, H ¼ f0; 0; Hg is the external magnetic field (always aligned along the OZ axis as per the previously discussed convention), H D is the demagnetizing field, created by the magnetization distribution MðrÞ, which explicitly depends on r and has a functional dependence on the whole MðrÞ vector field via the Maxwell equations, K 0 and KðrÞ denote the constant global and small spatially fluctuating local anisotropy parameters, s and dðrÞ are the corresponding anisotropy axis directions (jdj ¼ jsj ¼ 1), and 0 is the permeability of a vacuum. The anisotropy is of the easy-axis or easy-plane type, depending on whether K 0 > 0 or K 0 < 0. We also assume that the saturation magnetization is weakly fluctuating: where the magnitude of I M ðrÞ ( 1 is a small quantity with a zero spatial average hI M ðrÞi ¼ 0, so that M 0 ¼ hM S ðrÞi is the average saturation magnetization of the magnet. The quality factor of the fluctuating anisotropy ( 1 is assumed to be small and of the same order as I M with its average being zero: hI K ðrÞi ¼ 0. To cover several systems of magnetic units, the magnetic induction is defined here as (Aharoni, 1996). The model (4) with K 0 ¼ 0 has already been studied by Honecker & Michels (2013), by Metlov & Michels (2015) and in many other papers under a wide spectrum of assumptions about the random material parameter fluctuations (random defects). Global uniaxial anisotropy K 0 6 ¼ 0, which is new in (4) in the context of magnetic SANS, can have several physical origins. It may appear due to heat treatment in a magnetic field (magnetic annealing). It can be the result of an applied uniaxial mechanical stress or annealing of a material that is simultaneously subjected to an applied stress (stress annealing). Or, it can be imprinted into the material by severe plastic deformation (e.g. by high-pressure torsion).
The equilibrium magnetization vector distribution, minimizing (4) under the constraint jMj ¼ M S , is a solution to Brown's equations (Brown, 1963;Aharoni, 1996): where the cross denotes the vector product. The effective magnetic field H eff ðrÞ is defined as the functional derivative of the ferromagnet's energy-density functional EðfMðrÞgÞ over the magnetization vector field: Linearity of the variational derivative allows one to compute the contribution of each term in (4) separately. They have been published earlier (Metlov & Michels, 2015). Let us give here only the expression for the effective field associated with the magnetic anisotropy [the last two terms in (4)]: is the global anisotropy quality factor. Note that expression (8) is only valid up to the first order in small quantities I K and I M ; higher-order terms of this expansion are ignored. The contribution / 2I M in the last term appears due to the series expansion of M 2 S ðrÞ in the denominator of the corresponding term in (4). In this work, we will limit ourselves only to the research papers second-order theory (first-order solution of the micromagnetic problem, second-order contribution to the magnetic SANS cross sections) in the amplitude of small material parameter fluctuations. Consequently, the exchange length L E can also be considered as a nonfluctuating constant L 2 its small fluctuations only contribute to higher orders (Metlov & Michels, 2015).
When the material parameter fluctuations are small (I M ; I K ( 1), the approximate solution of the micromagnetic problem (4) can be obtained as a perturbation on top of the ground state that forms in the absence of fluctuations (I M ¼ I K ¼ 0). In our case, the ground state always corresponds to the uniform magnetization state, for which the first and the third terms in (4) are zero because they are proportional to the spatial derivatives of the magnetization. Minimization of the remaining second (Zeeman energy) and last (global anisotropy energy) terms defines the orientation of the uniform ground-state magnetization. In the case of K 0 ¼ 0 the magnetization is strictly parallel to the external field H. The global uniaxial anisotropy will rotate the magnetization away from the direction of the field H by the angle towards the direction of the anisotropy axis s, which can be parametrized via spherical angles s and s . The corresponding equation for the angle between the vectors M and H follows from (6): where the dimensionless magnetic field equals h ¼ H=ð B M 0 Þ. It is identical to the one studied by Stoner & Wohlfarth (1948) and in the general case may exhibit quite a nontrivial hysteresis. In the case when the anisotropy axis is oriented perpendicularly to the external magnetic field ( s ¼ =2) the lowestenergy solution of (10) can be written as where H K ¼ 2K 0 =ð 0 M 0 Þ is the critical field at which saturation (of the coherent rotation of magnetization) occurs and the average magnetization becomes strictly parallel to the applied field H. Another interesting case is when the anisotropy axis is directed along H ( s ¼ 0). In this case, a square hysteresis loop with a width of H K will be observed. Uniform rotation of magnetization by itself is irrelevant to the SANS experiment as it produces no scattering of neutrons, but its interplay with random material parameter fluctuations, as we will see later, makes a substantial impact on the magnetic SANS cross sections at q 6 ¼ 0.
Because the ground state of the considered magnet is always uniform, we introduce another (primed) coordinate system, in which the OZ 0 axis is oriented along the groundstate magnetization vector. Specifically, the original and the primed coordinate systems are connected by a rotation in the plane containing the magnetic field vector H and the direction of the global anisotropy s. Denoting the corresponding angle as , this rotation can be described by the matrixR RðÞ, which relates the coordinates of vectors in the original and primed coordinate systems, e.g. s 0 ¼R Rs. By construction, among the magnetization vector components in the rotated coordinate are of the first order in I M and I K and the last one can be found from the constraint jMj ¼ M S as M 0 Z ' M 0 þ M 0 I M up to the terms of the first order in I M .
Because the SANS cross sections (1) and (2) are expressed in terms of the Fourier components of the magnetization e M M, it is convenient to solve Brown's equations directly in Fourier space. This also makes the solution of the Maxwell equations for the demagnetizing field easier, but turns products of realspace functions in the effective field and in Brown's equations into convolutions. They are denoted here as The algebra of convolutions is commutative, distributive and associative with respect to multiplication by a constant. It also has an identity element such that È P ¼ P.
Using this notation we can express the Fourier image of the effective field (7) including all the terms from the energy (4) as where the convolution of two vectors ðd 0 È e M M 0 Þ is understood as their scalar product with multiplications replaced by convolutions. We also assume here, as did Metlov & Michels (2015), that e M M tends to 0 fast enough as q ! 0 that the third term in (13) tends to 0 in this limit. Expression (13) is only correct up to the first order in I M ; I K ( 1. To this order it is, basically, the same as the Fourier image of the effective field in the work of Metlov & Michels (2015) apart from the addition of the last term and the fact that it is expressed in the primed coordinate system, whose OZ 0 axis may rotate away from the direction of the external field H.
Of the three Brown equations for each vectorial component of (6) only two are independent. They also contain convolutions in Fourier space: Substituting the components of the effective field and the Taylor expansion of the magnetization M 0 ¼ fM 0ð1Þ X ; M 0ð1Þ Y ; M 0 þ M 0 I M g (where the quantities M 0ð1Þ X , M 0ð1Þ Y are of the same order as I M , I K ), Brown's equations become Taylor series themselves.
To zero order we recover equation (10), which has already been analyzed. Collecting the first-order terms leads to the following system of equations for the components of the normalized magnetization vector m 0 ¼ M 0 =M 0 : research papers components of the q vector are fx q 0 ; y q 0 ; z q 0 g ¼ q 0 =q. This linear system of equations can be easily solved: where W ¼ s 0 X y q 0 À s 0 Y x q 0 . These two expressions are only different with respect to exchanging X $ Y everywhere (note that W changes sign under such a transformation).
In the absence of global anisotropy (Q ¼ 0) the primed coordinate system is identical to the original one. Then 0 q 2 and the resulting expressions for e m m 0ð1Þ X and e m m 0ð1Þ Y coincide with the ones obtained by Metlov & Michels (2015), which in turn have several limiting cases in the earlier micromagnetic SANS theory and in the approach-to-saturation theory.

Magnetic SANS cross sections in the low-field limit
In this section, we analyze the perpendicular scattering cross sections in the regime of low magnetic fields (H < H K ). For brevity and for practical relevance (as illustrated in Section 6) only the particular case when the anisotropy is directed along the OY axis is considered.
When the anisotropy is directed along OY the equilibrium magnetization deviates from the field direction by rotating in the YOZ plane around the OX axis. This means that h 0 X ¼ 0 and s 0 X ¼ 0. Assuming that the anisotropy and saturation magnetization inhomogeneity functions are related via a scalar factor , i.e. I K ¼ I M , and setting I M ¼ I, noting that in the perpendicular SANS geometry x q 0 ¼ x q ¼ 0, substituting (18) and (19) into (1), and performing an averaging procedure over the defect realizations and the representative volume orientations described by Metlov & Michels (2015), one can obtain the following expression for the macroscopic magnetic spin-misalignment SANS cross section: h e I I 2 i is the average of the squared Fourier image of the inhomogeneity function e I I, which for the Gaussian randomly placed inhomogeneities were calculated by Metlov & Michels (2015).
Note that the cross section is proportional to h e I I 2 i, because linear terms in I K and I M average to zero. We emphasize that the presence of the anisotropy in the OY direction rotates the SANS cross section as a whole and (20) gives the scattering cross section in the rotated coordinate system as a function of q 0 ¼ qf0; y q 0 ; z q 0 g. In the perpendicular scattering geometry [see Fig. 1(a)] at low fields (h < Q), this rotation by the angle (11) in the detector plane can be described using the following matrix:R and In particular this means that in the absence of an applied field (h ¼ 0) the s 0 vector is oriented along the OZ 0 axis of the primed coordinate system and the average magnetization is directed along the easy axis. By contrast, in a saturating field (h ¼ Q), the field h 0 is along OZ 0 , implying that the magnetization is parallel to the applied field. One can obtain dAE ? M =dðq Y ; q Z Þ (note the absence of primes) from dAE ?0 M =dðq 0 Y ; q 0 Z Þ by substituting Formally, expression (20), prior to substitution of the rotation matrix, is valid in both the high-and the low-field limits with an appropriate selection of the primed coordinate system. In the high-field limit, it can be further simplified as described in the next section.

Magnetic SANS cross sections in the high-field limit
In the saturation regime (H > H K ), the macroscopic mean magnetization is directed along the external magnetic field. This means that there is no need to rotate the coordinate system and the expressions for the cross sections become simpler. (20), and removing all the primes one obtains the following expression:

Perpendicular SANS geometry
In the limiting case of Q ¼ 0 and after introduction of the polar angle for the scattering vector in the detector plane, q ? ¼ qf0; sin ; cos g, it coincides with equation (49) of Metlov & Michels (2015). It is convenient (Honecker & Michels, 2013) to decompose the spin-misalignment SANS cross section into the contributions related to anisotropy and saturation magnetization fluctuations, where dAE res =d is the (nuclear and magnetic) residual SANS cross section that is measured at complete saturation and is assumed to be magnetic field independent. The function S H ðqÞ ¼ Vb 2 H M 2 0 h e I I 2 i2 2 =15 is called the anisotropy-field scattering function, whereas S M ðqÞ ¼ Vb 2 H M 2 0 h e I I 2 i is the scattering function of the longitudinal magnetization. Because of the averaging over the scattering volume orientation these functions depend only on the magnitude of the scattering vector. The micromagnetic response functions R ?;Y H and R ?;Y M , whose superscript here marks the anisotropy axis orientation, have the following form: There is no new angular dependence here, but for Q 6 ¼ 0 the relative strengths of the isotropic halo and the scatteringangle-dependent terms are modified. (24)] decomposition can also be introduced in this case with the following response functions:

s k OX. A similar [to
The global anisotropy does not affect R ?;X M , only R ?;X H . 5.1.3. s k OZ. When the anisotropy axis is oriented along the direction of the magnetic field, the magnetic SANS cross section is given by Its only difference from the globally isotropic case (Metlov & Michels, 2015) is the replacement of h q by h 0 q , which means that the anisotropy only modifies (increases or decreases, depending on the sign of Q) the effective external field strength.

Parallel SANS geometry
Usually the parallel magnetic SANS cross section is isotropic in the detector plane. However, the presence of a global anisotropy breaks this property. The cross section can be computed by substituting (18) and (19) into (2) and averaging over the defect realizations and the representative volume orientations. This yields where the response functions R k;X=Y H are no longer isotropic and depend on the polar angle of the scattering vector q k ¼ qfcos ; sin ; 0g. They are given by the following expressions: where the first equation refers to the case s k OX and the second expression is for the s k OY case. Generally, the SANS cross section is compressed (Q > 0) or expanded (Q < 0) along the anisotropy axis. When the anisotropy axis coincides with the direction of the external magnetic field (s k OZ), the parallel spinmisalignment SANS cross section remains fully isotropic in the detector plane, but the denominator now contains h 0 6. Analysis of experimental data 6.1. Field-annealed Vitroperm One way to induce a global anisotropy into a magnetic nanocrystalline material is to anneal it in a magnetic field. Michels et al. (2003) and Grob et al. (2004) report magnetic SANS data of a field-annealed Vitroperm alloy, which is a twophase iron-based nanocrystalline soft magnetic material. The azimuthally averaged perpendicular spin-misalignment cross section data are shown in Fig. 2 as a function of the magnitude of the scattering vector q for a number of different values of the applied magnetic field H. The direction of the anisotropy axis is shown in the inset of Fig. 2. From the Vitroperm hysteresis loop, shown in Fig. 1 of Grob et al. (2004), we can research papers determine that magnetic saturation by coherent rotation occurs at about 40 mT. This means that among the available data the four fields of 0.9, 11, 14, and 24 mT correspond to the low-field regime, considered in Section 4.
The data for the remaining three field values correspond to the high-field regime, but, unfortunately, are not sufficient to test the corresponding high-field cross section expression. This is because the procedure of fitting the data by (25) involves linear regression to determine the three unknowns (the residual cross section and the scattering functions S H and S M ) for each value of q, based on (three) cross section data points, corresponding to each field value. It is always possible to draw a plane through any three points in 3D space, which means that the fitting would not help to validate the high-field cross section model. This case is tested separately in Section 6.2 on a different data set with many more high-field data.
In terms of the response functions (25) the spin-misalignment SANS cross section in the low-field regime in the primed coordinate system (20) can be represented by To interpret the azimuthally averaged cross section data, the response functions need to be averaged over in q ? as well by computing h. . .i ¼ 1=ð2Þ R 2 0 ð. . .Þ d. The rotation of the cross section as a whole in the plane of the detector is insignificant for the azimuthal averaging, so that we can integrate directly in the primed coordinate system by setting y q 0 ¼ sin and z q 0 ¼ cos . Substituting the values of h 0 and s 0 from (22) and performing the integration we get where ¼ L 0 q and g ¼ Àh 2 þ QðQ þ 2 Þ. Taking the known value of the exchange stiffness A = C/2 = 1 Â 10 À11 J m À1 (or L 0 = 4.2 nm) and the saturation magnetization 0 M 0 = 1.2 T of Vitroperm, only the parameter Q in the above functions is unknown. To determine Q, we have performed a linear regression of the cross section data for the four chosen values of the field corresponding to the low-field regime. At each q, the data were fitted by expression (25), obtaining dAE res =dðqÞ, S H ðqÞ and S M ðqÞ, and then the total least-squares error of this fit was numerically minimized to obtain the value of Q. This procedure yields Q ¼ Q VP ¼ 0:0277 AE 0:0002 and the fitted curves are displayed by the solid lines in Fig. 2. Note that the field dependence of the fitted curves is solely determined by the explicit field dependence of the average response functions hR ?;Y H i and hR ?;Y M i . The linear regression parameters dAE res =d, S H and S M are field independent. The lower-left inset in Fig. 2 shows the total error of the fit around the optimum value of Q ¼ Q VP .
To verify the self-consistency of the fit, we have repeated the procedure fixing the value of Q ¼ Q VP and treating L 0 as an adjustable parameter. The minimum of the total leastsquares error (shown as an inset in Fig. 2) corresponds to L 0 = 4.039 AE 0.002 nm. The errors were computed using a Monte Carlo procedure by adding a random AE5 cm À1 sr À1 contribution to the measured dAE=d values and computing the standard deviation of the resulting Q and L 0 across many realizations of this random process.
Using the obtained value of Q, we can compute the value of the uniaxial anisotropy as K 0 = 0 M 2 0 Q=2 ffi 15 900 J m À3 , which is much larger than the 10 J m À3 estimated by Grob et al. (2004) from an analysis of the domain-wall width. However, when the external field is perpendicular to the easy axis, the saturation field is of the order of the anisotropy field. Using the obtained value for K 0 yields 0 H K ¼ 2K 0 =M 0 = 33 mT for the anisotropy field, which is in reasonable agreement with the saturation magnetic field > $ 10-20 mT measured by Grob et al. (2004). Such large uniaxial anisotropy values were previously also reported by Herzer et al. (2011) for Vitroperm samples annealed under tensile stress.

High-pressure-torsion nickel
Due to its inherent axial symmetry, high-pressure torsion (HPT) is another way to induce a uniaxial anisotropy in an originally isotropic sample. Oba et al. (2021) report magnetic SANS cross section data (shown in Fig. 3) on an Ni sample subjected to such a treatment.
In the HPT-Ni experiment, the natural anisotropy axis of the sample (the axis of torsion) was aligned with the neutron beam during the SANS experiment and an isotropic parallel SANS cross section was observed. It is well known that this isotropy indicates that the sample itself is macroscopically Azimuthally averaged cross section data ðdAE ? =dÞ M of Vitroperm at selected applied magnetic fields (log-log scale); the experimental data are those from Michels et al. (2003). Solid lines show the fit by (25), which is also valid in the low-field limit, using the response functions (33). The insets show the total mean-square deviation error for the Q and the L 0 fits described in the text, as well as the relative orientation of the anisotropy axis and the applied magnetic field.
isotropic. Yet, as we can see from Section 5.1, despite the fact that there are no new angular terms when the neutron beam is parallel to the anisotropy axis, the cross section is still slightly modified. These changes can be extracted from the experimental data together with the corresponding anisotropy constant (quality factor Q) value.
We have fitted the spin-misalignment perpendicular SANS cross section of HPT-Ni (at all q and H values besides 0 H = 0.1 T as will be explained later) from Oba et al. (2021) with the response functions from (27). The sin 2 term is introduced and discussed by Oba et al. (2021) and is beyond the scope of the present work. Also, Oba et al. (2021) assume S H ðqÞ and S M ðqÞ depend on the entire q vector, meaning that each angle can have its own fitted values of the scattering functions. Here, we assume (as follows from the directional averaging procedure) that these two scattering functions depend only on the magnitude of the scattering vector and the angular dependence of the cross section (34) is entirely contained in the response functions. Even with such a restriction (and consequently much less fitting freedom) we were able to obtain a fit of comparable quality [see the solid lines in Fig. 3 and in Fig. 8 of Oba et al. (2021)] by minimizing the total error with respect to the Q value in the response functions. The fit with Q ¼ 0 (dashed lines in Fig. 3) is significantly worse, especially at the low values of q, where the cross section values are the largest. There is a systematic leaning of the cross section curves towards smaller , especially pronounced for smaller q, which we attribute to a possible slight misalignment of the anisotropy axis from the direction of the neutron beam. In the analysis of Oba et al. (2021) this systematic leaning was absorbed into the angular dependence of the scattering functions S H ðqÞ and S M ðqÞ. The best-fit value for the anisotropy quality factor is Q = Q HPT-Ni = 0.18. For this reason, we had to omit the 0.1 T data from the fitting procedure, as (unlike the rest of the data) they do not fall into the high-field regime for which the response functions (27) were derived. The anisotropy value turns out to be rather large, K 0 ¼ 0 M 2 0 Q=2 ffi 2.6 Â 10 4 J m À3 (using M 0 = 482 kA m À1 ). It was silently absorbed into the S H ðqÞ and S M ðqÞ response functions in the work of Oba et al. (2021) but can be revealed using the present more sophisticated theory; and it agrees very well with the value of 1.8 Â 10 4 J m À3 estimated by Bersweiler et al. (2021) using a correlation-function analysis on the same specimen.

Summary and conclusions
The existing micromagnetic SANS theory for spatially inhomogeneous ferromagnets is extended here by including the effect of a nonzero average uniaxial anisotropy with quality factor Q. The anisotropy leads to a deviation of the average magnetization from the external magnetic field direction accompanied by the Stoner-Wohlfarth rotational hysteresis. The macroscopically averaged (over the orientation and realizations of the random material defects in the scattering volume) SANS cross sections are computed analytically, on the basis of the presented solution of Brown's equations of micromagnetics in the small-misalignment approximation. In view of their simplicity and practical significance, we have compiled the cross section expressions for several special cases, where the anisotropy axis is either perpendicular or parallel to the magnetic field direction in both the parallel and the perpendicular SANS geometry. It follows from Stoner-Wohlfarth theory that for the mutually perpendicular anisotropy axis and the magnetic field direction there is a critical magnetic field at which saturation of the coherent rotation of magnetization occurs. We have analyzed the SANS cross sections both below and above this rotational magnetic saturation and computed the micromagnetic SANS response functions for the latter case. Some of these expressions exhibit an additional angular dependency of the scattering cross section on the scattering vector orientation, compared with their previously known Q = 0 limit. The present theory fits well the azimuthally averaged SANS cross section of fieldannealed Vitroperm alloy and the angular dependence of the spin-misalignment SANS cross section of high-pressuretorsioned Ni, allowing determination of the respective global anisotropy quality factors.